Vacuum Rabi oscillation in nonzero-temperature open cavity 



(N 

O 

(N 



Oh; 



> 



O 



X 



Patrycja Stefanska,^ Marcin Wilczewski,^ and Marek Czachor^ 
^ Zespdl Fizyki Atomowej — Katedra Fizyki Atomowej i Luminescencji, Politechnika Gdanska, 80-233 Gdansk, Poland 
^Katedra Fizyki Teoretycznej i Informatyki Kwantowej, Politechnika Gdanska, 80-233 Gdansk, Poland 

Published as: Open Syst. Inf. Dyn. 18 (2011) 363-373 
DOI: 10.1142/S123016121100025X 

Comparison of tlieory of Rabi oscillations with experiment [M. Wilczewski and M. Czachor, Phys. 
Rev. A 79, 033836 (2009)] suggests that cavity lifetime parameters obtained in measurements with 
many photons may be much smaller than those applicable to almost vacuum states of light. In this 
context we show that the conclusion remains unchanged even if one takes a more realistic description 
of the initial state of light in cavity. 

PACS numbers: 42.50.Lc, 42.50.Dv, 32.80.Ee, 32.80.Qk 



I. INTRODUCTION 

The 1996 Brune et al. experiment on vacuum Rabi os- 
cillation [1, 2] was analyzed in Q by means of several 
alternative models of atom-reservoir interaction. The 
study was motivated by difficulties with fitting the data 
by theoretical curves, a problem addressed earlier by var- 
ious authors 0-01 • Agreement with experimental Rabi 
oscillation data was then obtained but for the price of 
a cavity quality factor that was 500 times bigger than 
the one reported in [1]. A part of open questions thus 
remained. 

The solutions of master equations discussed in Q were 
easier to find than in standard approaches because the 
formalism was based on jump operators generating tran- 
sitions between the dressed states of the atom-field sys- 
tem. Such a construction is more consistent with the 
general theory of open systems 0,13 than the popular ap- 
proach based on jumps between the atomic (hence bare) 
states [3], and is mathematically simpler. In the context 
of quantum optics it appeared only relatively recently in 
[To']. Another reason why it was possible to find exact 
solutions was that the atom-field system was assumed to 
start from the initial photonic vacuum state. In terms of 
dressed states the initial state belonged to the subspace 
of the first dressed-state doublet. 

The latter assumption was not very realistic. Light in 
the cavity was initially in thermal state at T = 0.8K. 
The analysis from Q not only took into account emis- 
sions from the dressed states downward on the energy 
ladder, but also thermal excitations from vacuum to the 
two dressed states, as well as thermal long-wave fluctu- 
ations between the dressed states from the first doublet. 
But it did not take into account the presence of the re- 
maining bands of dressed states in the initial state. The 
first doublet appears with probability around 0.95, but 
the second doublet has probability higher than 0.045. It 
looks like downward transitions from the second doublet 
are processes of the same order as the upward thermal 
transitions at 0.8K from exact vacuum to the first dou- 
blet. So what remained unclear in Q was to what extent 
the fact that the initial state was thermal was influencing 
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FIG. 1: Probabilities of vacuum (n = 0) and 1-photon initial 
states (n = 1) effectively dominate the initial thermal prob- 
ability distribution at 0.8K. Contributions from n > 1 were 
not taken into account in 



the behavior of the system. In order to understand the 
issue one could iteratively solve equations involving more 
and more doublets and more transitions between them, 
and compare predictions with the data. If inclusion of a 
next doublet would not produce visible modifications of 
Rabi oscillations, it would be justified to conclude that 
truncation of the Hilbert space to a subspace spanned 
by a given number of dressed state has sufficiently well 
approximated the thermal state. 

The goal of the present paper is perform this test on the 
data from Jj] . We will see that inclusion of the next dou- 
blet of dressed states essentially complicates calculations, 
but does not really change agreement with experiment. 
We conclude that solution of the problem of the "wrong" 
cavity Q factor will not be achieved by taking more re- 
alistic initial states. The true physical mechanism must 
be therefore different. 



II. OPEN-CAVITY MODEL 

We employ the standard Jaynes-Cummings Hamilto- 
nian H ~ hfl in exact resonance. 
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FIG. 2: Energy levels and decay coefficients used in the gener- 
alization of the Scala model discussed in E± — hQ±, 



The initial state at T = 0.8 K is 



p(0) = po|e,0)(e,0|+pi|e,l)(e,l| + .. 



(2) 



where po = 0.952381, pi = 0.0453515, p2 = 0.00215959. 
Let us note that Y^'^=iPi = 1 ^ Po = 0.047619 is of the 
order of pi. The analysis performed in [3] assumed jumps 
between dressed states with transition coefficients typical 
of T = 0.8 K, but the initial state was approximated by 
p(0) — |e, 0)(e,0|. The solution given in employed 
transitions between the dressed states shown in Fig. 2, 
i.e. those involving |e,0) and the ground state: 



\^+) = ^(|.g,l> + |e,0)), 
l^^-> = ^(|g,l)-|e,0)). 
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The decay coefficients from Fig. 2 satisfied 
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e ra— 71 « 0.046632771, 



■72 ~ 0.046632872, 
e"W73 ~ 0.99999773, 
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Thermal excitations from the ground state to the first two 

dressed states involve proportionality factors e 
that are of the same order as pi. Since pi measures 
the probability of occurrence of |e, l)(e, 1| in the initial 
thermal mixture, it simultaneously determines probabil- 
ities of finding the next two dressed states. Accordingly, 
transitions determined by 7a and 7^ should be regarded 
as processes of the same order as downward transitions 
from the second dressed-state doublet: 



|f^2+) = -^(|<?,2) + |e,l)), 
1^2-) = -i=(|<?,2)-|e,l)). 



(9) 
(10) 



It is therefore justified to regard Fig. 3 as more realistic 
than Fig. 2. The open-cavity generalization of the Scala 
et al. model [l^l, corresponding to Fig. 3, reads 
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Expressing the initial condition p(0) = J2j '^jPj ^ 
ear combination of eigenvectors of £, Cpj — A^jPj, we 
find the solution p{t) — cje^'^pj. We need 25 eigen- 



vectors. Twenty of them are easy to find and determine 
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FIG. 3: E2+ = m2+,E2- = nn2-, e± = hQ±, Eo = mo. 



off-diagonal matrix elements (in the dressed-state basis), 



P6 = |^^2+)(^^2-|, 

P7 = |02+)(fi+|, 

Ay = -^(^2+ - ^^+) - 



(12) 



74 + 75 + 76 + 77 + 78 + 7e 
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i{n2- - 0+) - - 



71 + 73 + 74 + 75 + 76 



7c + 72 + 74 + 75 + 76 



7a + 76 + 74 + 75 + 76 



7o + 76 + 77 + 78 + 7e 



72 + 77 + 78 + 7c + 7e 



71 + 73 + 77 + 78 + 7e 
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71 + 72 + 73 + 7c 



71 + 73 + 7a + 76 



72 + 7a + 76 + 7c 



7i + 72 + 73 + 7c 



72 + 77 + 78 + 7c + 7e 



7c + 72 + 74 + 75 + 76 



7a + 76 + 74 + 75 + 76 



7a + 76 + 77 + 78 + 7e 
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72 + 7a + 76 + 7c 



The remaining five eigenvectors pj, j = 1,2,3,4,5 are 
related to the diagonal matrix elements of p, 



Pj = Xj\n2+){n2+\ + yj\n2-){n2-\ + zj\n+){ii+\ 
+vj\n-){n-\ + wj\no){i^o\ (33) 



The corresponding eigenvalue problem is equivalent to 
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The eigenvalues are 

Ai = 0, 

A2 -0.25 + (5 + \/^) , 

A3 = -0.25 fcj + 5 - V^V 



■where 



A4 = -0.25 (c + e + \/^) 

A5 = -0.25 (C + C - V^) 



(35a) 
(35b) 

(35c) 

(35d) 
(35e) 



K = 76(('5+w)2 -4(72(73 + 7a) (3g^^ 

+ (7a +7b)(73 ■4^7c) +71(72 +76 + 7c))), 

^ = (C-0'+4757e (36b) 
w = 71 +72 +73, C = 74 +75 +76, (36c) 

^ = 77 + 78 + 7e, (5 = 7a +7f) + 7c- (36d) 

Exphcit forms of the eigenvectors can be found in the 
Appendix. 

A reasonable approximation of the initial thermal state 
is given by 

p(0) = 0.95|e,0)(e,0| +0.05|e,l)(e,l| (37) 

= o.95(|r2+)(i7+| + |r2-)(r2-| 

+ 0.05(|r!2+)(f^2+| + |f^2-)(f^2-| 

- |l]2+)(f^2-|-|f^2-)(f^2+0. (38) 



After somewhat lengthy but simple calculations one finds 
that 



p{t) = 2{Bie'^''pi + B2e''^'p2 + B^e^'^'p^ + B^e'^^'p^ 
^^^'p,)+B,e^^'f 

„Al6t„ I D c^iot. 



B,e''''p5) + B^e'^^'pQ + ^136'^- Vi3 



-Bi6e^^i«Vi6 + ^i9e^^^«Vi9- 



(39) 



The coefficients are explicitly given in the Appendix. 
Probability pg{t) = Pgfi[t) + Pg,i{t) + Pg.2{t) of finding 
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the atom in its ground state reads finally 

Pg{t) = Bi{xi+yi + zi+vi + 2wi)e'^^* 

+ B2{X2 + y2 + Z2 + V2 + 2w2)e^ 

+ B3{x3 + y3 + Z3 + V3 + 2w3)e^^' 
+ B4{x4 + y4 + zi + vi + 2w4)e^''* 
+ B5{X5 + y5 + Z5 + V5 + 2w5)e^'-' 

74+7'5+'76+7'7 + "^S+7e + , — 

- 0.025e 2 * cos 2y2gt 

71 +7-2+7'3+7'c 



0.475e^ 



3-r73_Z7c+ 1— 

*cos2V2gt. (40) 



Let us stress that the above solution is found under the 
assumption that the atom-field coupling g is constant in 
time. We know, ho-wever, that the atom interacts with 
the mode whose spatial profile is Gaussian with width w. 
The atom is propagating through the cavity which makes, 
effectively, the coupling time-dependent. A method of 
taking this into account was discussed in detail in 
Assuming that the length of the cavity is d we obtain 
probability appropriate for comparison with experimen- 
tal data 

Pg(t) = Bi(a;i +yi + zi +wi -f 2u;i)e'^i* 

+ ^2(^2 + y2 + ^2 + W2 + 2?i;2)e'^"* 
+ B3{x3 + y3 + Z3 + 173 + 2?i;3)e^='* 
+ 54(2:4 + y4 + 2:4 + W4 + 2^4)6^"* 
-h 55(2:5 + ys + 2:5 + W5 + 2w5)e^«* 
- 0.025e 1 * cos (2%/2gV^-t) 



0.475e" 



71+72+73+71 



-*cos(2\/2g^^F^t). (41) 



Let us recall that the data shown in jll were plotted as 
a function of an effective time igff — V^^'^- Since it is 
more convenient for us to work in terms of t than igff 
(sf, we rescale the data to t. Fig. 4 shows vacuum Rabi 
oscillation predicted by our more realistic scenario (solid 
line) as compared to experimental data and predictions 
of the simplified model from ^3*1 (dotted). 



III. CONCLUSIONS 

The effort devoted to solving the more realistic case 
apparently did not pay: The new curve does not describe 
the data any better. The approximation made in Q turns 
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FIG. 4: The solid line represents (|41|l . The dotted curve is 
the prediction from 0]. The parameters are 71 = 72 = 74 = 
76 = 77 = 78 = 17.73 Hz, 7a = 75 = £71, 73 = 7c = 75 = 7e = 
0.07g, g = 477rl0^ Hz. t is the true time. 



out to be physically reasonable. This is in a sense good 
news since inclusion of yet higher dressed doublets would 
require solving algebraic equations of order higher than 5, 
and thus there would be practically no chance for finding 
exact solutions. 

On the other hand, this also means that refined the- 
oretical analysis of Rabi oscillation has not brought us 
any closer to understanding why probability of atomic 
exited state does not decay to zero as fast as expected 
on the grounds of cavity lifetime reported in The 
cavity quality seems to be better than that assumed by 
Brune et al. Perhaps, as suggested in Q, the key el- 
ement is in opening of the cavity and inclusion of the 
long-wave transitions within a given doublet of dressed 
states. The appropriate jump operators occur if one mod- 
ifies the cavity-reservoir coupling. The required interac- 
tion with the reservoir has stronger dependence on the 
number of photons inside of the cavity. Measurements of 
cavity lifetimes are typically performed with more cav- 
ity photons than in Rabi oscillation experiments. It is 
possible that the same cavity has a much longer lifetime 
if only a few photons are present. The problem requires 
further experimental studies. 



Appendix 

Coefficients introduced in previous section 

are: = B13 = -0.025 and Big = Big = -0.025, 

whereas for n = 1, 2, 3, 4, 5 we have: 

B„ = i[0.2375(A3n + ^4„)+0.0125(Ai„ + A2„)] 
X 

where 
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{i,j,i,™} = i 
{i,j,l,m}^k 



^ ^ SijkmXi Uj Zk Wm {I — 1, .-j 5) 

{i,j,k.m} = l 

{i.j,k,m}^l 



A^m — 



5 

E 

{i.j ,k,l}^m 



Sijki Xi Vj ZkVi (m = 1, .., 5) 



X = E ^ijklm Xt Uj Zk Vl Wm 

{ij,k,l^m} — l 

For eigenvalues (fc — 1,2,3,4,5) we define eigenvec- 
tors with coordinates {xk,yk, Zk,Vk,Wk} 

xi = 
2/1 = 

Zl = 727a + 7c 7a + 76 7c 

Vl = 737a + 71 7& + 737& 
wi = 7i72 + 7273 + 7i7c 

For fc = 2, 3 we have 



Xk = 








Vk = 








Zk = 3271 




K- 




Vk = 3271 


[av - 






Wk = 3271 


{aw 


+ (- 




4, 5 we have 








Xk = I67176 









Vk 



{i,k,l.m} = l 

{j,k^l,m}^i 



32717576 (cy + (-l)'=%/^d 
Zk = 327176 (c, + (-l)'=V0 4) 
Vk = 327176 (c„-f(-l)'^V^rf.„) 
Wk = 1287176 (c„ + (-l)'=v^d^) 
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where 



a-z = 275(7677 - 7478) ('^^ - S"^ - 272 ((5 + uj) 

+ 4(727a + 7c 7a + 767c)) +ij6iS + 0J - 2^) 

- 27573) (4767c + 2767c(w -6) - 47270(74 + 76) 

+ 74(15^ -wS + 2((5 + w)-47c(7a + 76))) 

- «;76(74('^ + 3w - 20 - 2(7274 + 7577 + 767c)) 
av = 475(7478 - 7677) (73 ('^ - <5) + 2(737c - 71 76)) 

- K76 (76 ((5 + ^-20 - 2(7374 + 7578 

- 76(72 + 7c))) + ((S + oj~ 2^76 - 27578) 
X (274(73(0; - (5) + 2(737c - 7i7b)) 

+ 76((7i + 73 - 7a - 76)^ - (72 + 7c)^ + 47l7a)) 

dw = 2((76((5 + w - 2^ - 27578)(7i74((5 + w) 

+ 7276(<5 - w) - 272(7174 + 7374 - 7276) 

- 271 7c(74 + 76)) + 275(7677 - 7478) 
X (2(7172 +7372 +7i7c) -71 (<5 + w)) 
+ «;76(7i74 +7276)) 

bz = 476((7i+73)(74(^-w)+ 7274 + 7577) 

- 7a(7l74 + 7276) + (w - C + 7c)767c 

- 7c (7374 +7578)) 

K = 476(7374(w-^ + 7c)- 76(7174 + 7276) 

+ (7578 + 76(C - 72 - 7c))(72 + 7c) 

- 73(7577 + 767c)) 

bw = 476 (7276 (<^ - - 72(7174 + 7374 - 7276 

+ 7578) - 7i(74(? - (5 - w + 7c) - 7577 - 767c)) 

cx = vic'-e)+o{e+r,-e 

+ C(4e - 4(5 - 4a; + 50) 

cy = e{2{S + uj)-3{^ + 0)-v{^ + 



Cz = 472((C - 0(7476 - 767a) + 275(787a - 7776)) 

+ (^ + C-2(5 + 27c)((27577-74(C-0) 

X (^ + C - 272 - 27c) + 27c(76(C - - 27578)) 

- 61(74 (^ + 3C - 2(5) - 2(7274 + 7577 + 767c)) 

Cv = -6'(76(-2((5 + w- 72- 7c) +3C + -2(7374 

+ 7578)) + (-2(5 + C + e + 27c)(273((C - 

X (74 + 76) + 275(7e-0)-(C + ?-27i) 

X (76(C - - 27578)) - 47i((C - 

X (7476 - 767a) + 275(787a- 7776)) 

Cw = 74 + 7276) - (7172 +71 7c + 7273) ((C - 
X (74 + 76) - 27578) + 7577(272(71 + 73) 
7i(C + C - 27c)) - 727578(C + 

dy = 2((5 + w-C-0(C + C)-r?-e 

dz = 4(7577(^ + C - (5 - 72) + C(74(<5 - C + 72) + 767c) 

- 75(787c + 747e) - (727a + 7c(7a + 76))(74 + 76)) 

dv = 4(-75(7377 - 7278 + 767e) - (74 + 76)(737a 

+ 76(71 + 73)) + 7578(C + C - (5 - w + 7c) 

+ C(7374 + 76(i5 + w - C - 72 - 7c))) 

dw = 7276 (C - w) - 72(7374 - 7276 + 7578) 

+ 7i74(C - 72 - 7c) - 71(7577 + 767c) 

V = 72(473 -2(^ + C- 27a)) 

+ (^ + C-273-27c)(e + C-2(5 + 27c) 

+ 271(272 -(^ + C- 2(5 + 27„)) 
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